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Abstract 



KA , We consider a Dirac field in 2 + 1 Euclidean dimensions, in the pres- 

5— 1 J ence of a linear domain wall defect in its mass, and a constant elec- 

tromagnetic field. We evaluate the exact fermionic determinant for 
the situation where the defect is assumed to be rectilinear, static, and 
the gauge field is minimally coupled to the fermions. We discuss the 
dependence of the result on the (unique) independent geometrical pa- 
rameter of this system, namely, the relative orientation of the wall 
and the direction of the external field. We apply the result for the 
determinant to the evaluation of the vacuum energy. 
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1 Introduction 

In general, the fermionic action S is (or may be transformed into) a quadratic 
expression in the Grassmann fields, which is in turn determined by a Dirac 
operator T>. This Dirac operator, except for some trivial cases, carries a 
dependence on some parameters and fields. At some stage in the resolution 
of a dynamical problem, it may be useful to regard those fields as 'external', 
either because they don't have a proper dynamics, or because they have not 
yet been quantized. 

The fermionic determinant is formally given by the product of the eigen- 
values of the Dirac operator. The explicit calculation of those eigenvalues is, 
however, impossible to achieve for an arbitrarily general situation, and one 
has to resort to some form of approximate expansion. Nevertheless, there 
are many highly symmetrical external field configurations where this explicit 
diagonalization can be successfully carried out. In particular, the constant 
F^ u case in several spacetime dimensions, allows for the determination of 
the exact eigenvalues, since the diagonalization of T> reduces in this case to 
finding the spectrum of a one dimensional harmonic oscillator |]J . Many in- 
teresting results have been found in this area, and they have application, for 
example, to the determination of effective actions ||. 

In most of these cases, the method of resolution amounts to expressing 
first the Dirac operator D as a quadratic form in the operators x^, p^ [] 
(as in the Fock-Schwinger method |1|]). The diagonalization procedure is 
then equivalent to defining a transformation from x^ and p^ to a new set of 
canonical operators, and thus may in general be expressed as a Bogoliubov 
transformation in some suitable creation and annihilation operators. The 
original problem becomes then an algebraic one, with a complexity which 
depends of course on the kind of situation considered. 

In this article we will evaluate the fermionic determinant corresponding 
to a Dirac field in 2 + 1 dimensions, coupled to an external uniform electro- 
magnetic field F^u, and with a mass term which is a linear functional of the 
coordinates. This kind of mass term is a simple example of a situation where 
a domain wall like defect (of constant slope) is present, in this case in the 
parity-breaking mass term. We note that, if the electromagnetic field were 
parallel to the defect, the situation would fall into the well-known Callan- 



1 We work with Euclidean coordinates: x^ = (xq,Xi,X2), where xq denotes the Eu- 
clidean time. 



Harvey phenomenon [[|, by which a chiral fermionic zero mode is induced 
on the defect. We will, however, consider situations where the field has a 
different orientation. 

We will show that the modulus of the determinant can be found exactly 
for these configurations, what represents a non-trivial generalization of the 
well-known situation corresponding to a uniform electromagnetic field and a 
constant mass term. 

The Euclidean action S for this system is 

S = f d 3 x^V^ , (1) 

where T> =$ + i A + M, and we have absorbed the electric charge into 
the definition of the gauge field A. The Hermitian 7-matrices are in the 
irreducible 2x2 representation of the Dirac algebra: 

{ 7/i ,7„} = 28^. (2) 

The domain wall, defined as the region where the mass changes sign, will 
in our case have only two defining properties: its location and its slope. The 
latter determines the localization (or not) of the defect, and is quantitatively 
measured by the normal derivative of the mass along the curve of the defect. 
The situation we consider in this paper may be considered as an approximate 
treatment of localized defects; a general discussion on the localization of 
modes for this system may be found in [|J . 

The Dirac operator T> appearing in ([!]) is not the more suitable one to 
deal with the eigenvalue problem, since it is not Hermitian. It is convenient 
to consider, instead, a related operator 7i which is always Hermitian, has a 
complete set of eigenstates, and its eigenvalues are the squares of the Dirac 
operator's eigenvalues. This operator TC is of course 

n = v ] v , (3) 

which is indeed the operator defining the positive Hermitian part of T> in its 
polar decomposition. We can see from the structure of H, that for a mass 
linear in the coordinates, and with a constant electromagnetic potential, it 
will be a quadratic expression in x^ and p M . The mass will be regarded as a 
linear functional of the coordinates. 

This kind of configuration for the mass and the electromagnetic field 
can be characterized by two vectors in Euclidean spacetime. For the mass, 



vanishing along a line that passes through the origin (what is always possible 
by a proper choice of coordinates), can be defined in terms of a vector rj^, 

M = v^v , (4) 

so that rjn points in the direction normal to the plane that defines the defect 
hyperplane. Regarding the electromagnetic potential A^, a linear function 
of the coordinates suffices to generate a constant electromagnetic field -F M „, 
since the gauge potential can be written in a symmetric gauge as follows: 

A/j, = — — r^yXy . (5) 

We have found it convenient to use, rather than the antisymmetric tensor, 
F^, its dual F M , 

f \i = ~^€^vp-t 'vpi \y) 

where e M ^ p is the totally antisymmetric tensor in 3 dimensions. This puts 
both the defect and the gauge field on an equal footing, simplifying the study 
of the dependence of the fermionic determinant on the geometric invariants 
that may be built out of those external fields. 

Writing A p in terms of the dual of F^ u , A p = ^e^ up F p x u , one can present 
the Dirac operator parametrized by two constant vectors, thus summarizing 
all its dependence on the external parameters: 

V = ^ ■ d — -'j ■ x x F + t] ■ x , (7) 

where we have used the notation 

(a x b)^ = t^ vp a v b p , a-b = a^b^ . (8) 

The relative orientation between rj^ and F p covers all the freedom to 
describe the different geometrical configurations allowed within this model 
for the Dirac operator. Different particular cases will describe rather differ- 
ent situations, both from the physical and the mathematical points of view. 
The simplest situation corresponds to rj parallel to F, as in the Callan and 
Harvey mechanism, since the gauge field is then entirely contained in the 
defect worldsheet. Then the determinant of the three dimensional system 
immediately factorizes into the product of an infinite number of (decoupled) 



determinants in 2 + 1 dimensions |3| . Our presentation begins, in section |2|, 
with the more interesting case F ■ r\ — 0. Namely, the electromagnetic field 
is normal to the defect hypersurface. The interest in this case stems from 
the fact that it takes into account the competition between the localizing 
effect due to the defect and the effect of the electromagnetic field, which will 
tend to 'drag' any charge distribution in a direction normal to the defect. 
An application to the calculation of the vacuum energy, for this case, is pre- 
sented in section ^|. Finally, we conclude by deriving the eigenvalues of the 
equivalent diagonal operator for the general case in section |j. 

2 Fermionic determinant for F ■ r\ = 

In this situation, F^ and rj^ define two orthogonal directions under the scalar 
Euclidean product defined in (§). Thus, in terms of these two vector fields, 
we may construct two families of curves in the Euclidean 3-dimensional space 
that may be used to define coordinates, varying along their integral lines. A 
third family of curves, orthogonal to the two previous ones, can then be 
uniquely defined by selecting a right-handed orientation. We define t^ to be 
the unit vector in the direction of F^. It is obvious from equation (^) that in 
the particular case we are dealing with, the operator T> does not depend on 
the coordinate x t (which parametrizes the integral curves of t M ) . We see that 
it is then sufficient to diagonalize the part of TC which depends on the other 
two coordinates. Indeed, using equation (0), we can express TC in terms of 
the external parameters, 

TC = -d 2 + -[x 2 F 2 -(x-F) 2 ] + (x-t 1 ) 2 - i(x x d) ■ F + 7 • (F - 77) , (9) 
which, in the coordinates defined by the unit vectors 

V - F 

fj = — '—, t = ^=-, b = fjAt, (10) 

WvW \\F\\' 



becomes, 



p2 p2 

5 2 x2 a2,/ r 1 2\ 2 , r 2 



U = -d 2 -d 2 -d 2 b + ( —+r ] 2 )x 2 ! + —xi-tF(x b d v -x v d b )+^-(F-r ] ). (11) 

In this expression, 7i is invariant under translations in x t , thus the diagonal- 
ization of TC amounts to solving the reduced problem defined by a different 



Hamiltonian H acting on functions depending on the two coordinates x v and 
x b , 

U = -d 2 t + 2H, (12) 

where 

H = ~vl + 2^6 + Y X * + lH + 2" (X ^ ~ XvPb) + 2 1 '^~ V) - (13) 
In fll3|) we have introduced the constants 



_p2 p2 

rf = T + V 2 , ^ 2 = T - (14) 

To study the diagonalization of H, it is first convenient to define suitable 
creation and annihilation operators for each coordinate, 

1,1/2 . -1/2 ,. 

a j = 7/fH X 3+Mj Pj), 



t 1/1/2 ■ -1/2 



V2 



Pj), J = V,b (15) 



(no sum over j), since the expression for H becomes more symmetric, and its 
diagonalization can be done by a generalization of the usual procedure for a 
quadratic Hamiltonian. 

In terms of these operators, H contains both bilinear and constant terms: 

_ j. + Lr , UJ„ uj b .. j. i. 

H = u v a ] v a v + u b ala b + — r— U /— - J — ){ a n a b ~ %%) 



Ai y u b y a; ' 



'/ 



\\F\\ [ujZ /ZJJT. . + + . 

+ -[^ + w 6 + 7.(^-77)], (16) 

which has the important property of containing terms which are quadratic 
in the creation and annihilation operators, and not just the standard one 
with creation and annihilation operators mixed by a Hermitian matrix. This 
makes the diagonalization more cumbersome. Indeed, to diagonalize H, we 
first introduce some 2x2 matrices, defined by 

H = a)Ajiai + -ajB jt ai - -a]Bjiaj + -[u> v + u> b + 7 • ( F ~ V)] , (I 7 ) 



where A is Hermitian while B is real and symmetric 

A=(^ ~ m ) B=(° + * 
y iu u b J y £ (J 

with the constants it and £ defined by 
it = 



(18) 





(19) 



The diagonalization is non trivial due to the presence of the B matrix, which 
only vanishes for rj = 0. To diagonalize the bilinear term in equation (|T7|) , we 
shall need to define new creation and annihilation operators, which will be 
linear combinations of the ones defined in Eq. [15]. We first note that the Oj, a, 
operators have been defined in such a way that they verify the commutation 
relations: 

,t J 



[a], a]} = [cij,ai] = 0, 

[H, a\] 

[H, ai ] 



a]Aji + icijBji, 
—ajA* t +iajBji, 



(20) 



where A* is the complex conjugate of A. We require the new operators bj, b, 
to be also canonical and to diagonalize the bilinear form of (|TT| ) . Following , 
we require them to satisfy the commutation relations, 

[6}, bj] = [bj, k] = 0, 
[bj,b]] =8ji, 

[H,b\]=b}E h 

[HM] = -biE l (21) 

(no sum over I), We can write the new operators in terms of the old ones by 
using a general Bogoliubov transformation of the following kind: 



bl 



aftji 



ajVji, 



bl = - a ]v; + aj u; t , 



(22) 



where, to fulfill the commutation relations above, the coefficients U and V 
must verify the matrix equations 

U ] U - V ] V = I , (23) 

U l V - V*U = , (24) 

AU - iBV = UE , (25) 

and 

iBU + A*V = -VE , (26) 

E being the diagonal matrix whose elements are the eigenvalues of H and 
W , as usual, stands for the adjoint of U. Similar redefinitions, although for 
a different case, where used in and ||. 

The system of equations involving the coefficient matrices defines a gener- 
alized diagonalization problem, which involves matrices rather than vectors. 
The next step is to obtain E explicitly in terms of the parameters of the 
theory. We first note that not all the equations that define this eigensystem 



are independent. To determine E we may, for example, use Equation (]25 ) to 
express the coefficients Vji in terms of Uji, the parameters of the theory and 
the eigenvalues E L . This yields 

v = ( riiC/n + s u U 21 r 12 U 12 + s 12 U 22 \ ^^ 

V r 21 U 2i + s 2 iU u r 22 U 22 + s 22 U 12 J ' 

where the coefficients Tji and Sji are defined by 

u u \ / i(Ei-w b ) i(E 2 -uJ b ) \ 

_U _U 1 5 S = I {(E-L-LJr,) i(E 2 -UJ v ) I ' (^8) 

t t J \ t t J 

Introducing this result into Eq. (|26|) yields two equations, now involving U\\ 
and U 2 i, plus two others involving U\ 2 and U 22 , 

2ioj n u 



Un = U 111Fir 



and 



E\ + t 2 — E'lCJfe — U 2 + -ElCJ^ — LUrjLUb 

El + t 2 - E x u n - u 2 + E x u h - u n uj b 
2u b u 



2iu b u 

U12 — —U 22 - 



' El + t 2 - E 2 u v - u 2 + E 2 u b - u v uj b 

El + t 2 - E 2 u) b - u 2 + E 2 uj v - u v uj b 
—lU 22 . [60) 



Demanding equations fl2"5|) and fl3"U|) to be consistent, implies a set of 
constraints for the eigenvalues E\. After some algebra, those constraints may 
be written in terms of just one equation, that determines the possible values 
of Ei in terms of the external parameters. There are 4 solutions to this 
equation: besides the double eigenvalue E\ = 0, we find: 



E 2 = -E 3 = y^ + 3u 2 b . (31) 

This set of solutions contains the true eigenvalues; however, to discard the 
spurious ones, we should check whether they are consistent with (|23|) and 



(|24|), which guarantee the canonical commutation relations for the 'new' op- 
erators. In terms of the variables defined in Eq.([27|), (124]) can be recast in 
the form, 

pq(s 22 - sn) + q{r 22 - r 2 i) + p(r 12 - r n ) + (s 12 - s 2 i) = 0, (32) 

where p and q are defined by: 

E 2 + t 2 - EjUb - U 2 + EjOJn - uo„uo b 

p = -i— — ■ — , (33) 



and 



2uo v u 



2u v u 



" El + t 2 - E k UJ b - U 2 + E k LU v - LU v LU b 

where the labels j and k stand for the two different eigenvalues that can be 
chosen from the four possibilities in Eq.(^TJ). 

Eq.(|23"l) splits into three algebraic equations, one corresponding to the 
non-diagonal elements, 

p + q* - (r n + qs n y(r 12 p + s l2 ) - (r 21 q + S 21 )*(r 22 + ps 22 ) = 0, (35) 

plus two others for the diagonal terms, 

\U U \ 2 [1 + \q\ 2 - |r u + q Sll \ 2 - \r 21 q + s 21 \ 2 } = 1, (36) 

|^ 22 r[l + br-|ri 2 p + si 2 | 2 -|r 22 +ps 22 | 2 ] = 1. (37) 

Eqs. (^) and (^) do not contain the coefficients of the matrices U and V, 
and thus they are equations for the eigenvalues E\. In terms of the original 
parameters, we obtain for the eigenvalues of H, 



#i = 0, E 2 = Ju 2 + 3u 2 = ^F 2 + r ] 2 . (38) 



Note that Q55| ) and (j3"7|) only fix the moduli of C/n and f/ 22 ; the phase of 
these coefficients is not fixed by the basis choice, and on the other hand no 
physical magnitude will depend on it. 

We note that one of the eigenvalues vanishes, what implies \Un\ — > oo. 
This means that, in the direction corresponding to the operators b± and b[, 
there is no harmonic mode, but rather a free motion. Thus we define the 
corresponding conjugate X\ and pi variables, such that Ti in (|TT| ) becomes 



n = P 2 t+Pi + 2E 2 b\b 2 + to v + uo h + 7 • (F - rf) . (39) 

We have thus obtained the diagonal form for the operator 7i. From this ex- 
pression, we can obtain the modulus of the determinant of the Dirac operator 
V. 



3 Vacuum energy in a constant electromag- 
netic field 

In this section we evaluate the vacuum energy for the fermionic system stud- 
ied in the previous section (i.e., when F-r] = 0); we shall see that this physical 
magnitude is entirely determined by (f|). 

To that end we first consider the vacuum energy density for a fermionic 
system in the presence of an external electromagnetic field. The vacuum to 
vacuum transition probability amplitude is given by the S matrix expectation 
value between vacuum states, which depends on the external potential A, 

S (A) =< in\S\0 in > . (40) 

This object is usually normalized with respect to the transition amplitude in 
the absence of external fields, S (A): 

MM 2 - gP£ . (41) 



In terms of the normalized transition amplitude S' Q (A), one then defines a 
local function w(x) || 

\S' {A)\ 2 = exp[- I d n x w(x)} , (42) 

10 



where n denotes the spacetime dimension. In the Euclidean formulation we 
can interpret w (x) as half the vacuum energy density. To see this we consider 
the Euclidean evolution operator in the interaction representation U(/3,Po), 
which obeys the equations 0, 

lim [Trtf(£ -£)] = (0,m|C/(oo,-oo)|0,m) = S' (A) (43) 

p^oo Z Z 

and 

hm[-ilnTr[/(§,-§)]=£ , (44) 

/3^oo p Z Z 

where Eq is the vacuum energy. From these equations and the definition 
(|42| ) we see that the integral of w(x) over the spatial coordinates is half the 
vacuum energy. In order to obtain w(x), we first have to evaluate Sq(A). In 
the functional integral representation, we may write 

S (A) = \M\ 2 I D^D^e- / rf3 -^^*, (45) 

where the notation T>[A] indicates the Dirac operator dependence on the 
external field A. Since the integral is over Grassmannian variables, we write 

S (A) = \Af\ 2 det(V[A}), (46) 

where det stands for the determinant over both spinor and functional spaces. 
Inserting ( |Hf) into the definition flilD , we see that 

\S' (A)\ 2 = exp{Tr \nH[A] - Tr \nH[0]}. (47) 

Now, we use the Frullani's identity: 

In - = lim/ —( e - sb -e- sa ) (48) 

b e-o I s 



to rewrite equation ( f47|) in a particularly convenient integral representation: 
log(|^(A)| 2 ) = Tr lim / -( e - sH W _ e ~ sH ^), (49) 

e^O J e S 

whence we can obtain the vacuum energy density as: 

™( e -sH[A] _ e sH[0}^ (5Q) 

11 



W(X) 



where now the trace only affects the spinor space indices. 

We evaluate now the vacuum energy, the integral of w(x) over all the 
Euclidean space, for the previously described case. The operators TC[A] and 
7i[0] are respectively given by, 

H[A] = p 2 t +p\ + 2b\b 2 E 2 + u n + u b + 1 -(F-T 1 ) (51) 

and 

H[0] = p 2 + M 2 - 7 • rj . (52) 

We can write TC[0] in terms of some creation and annihilation operators c 
and c\ defined as in equation (|15|), with frequency u = \\r]\\, as follows: 



HO] =p 2 t +pl + |M| (2c f c + 1) - 7 • V ■ (53) 

Evaluating the trace, we see that: 

2 lim f — J/ x \ e -slP 2 t+P 2 i+2blb2+^+u> b ]( e sE 2 _ e ->&)\ x \ 
e->0j e s I 

_ e -^P?+^+2ct C +|M|]^s|M|_ e -s||r,||^|^| _ (^ 



W[X) 



The system we are considering is not translation invariant, and moreover, 
we may write the explicit coordinate dependence of w(x) by using the com- 
pleteness of the eigenstates of the operators b 2 b 2 and c^c for the respective 
Hilbert spaces. Denoting by (x 2 |m 2 ) and (xjra,,) the respective eigenstates 
of the operators b 2 b 2 and c*c, which are as usual labeled by a non-negative 
integer m = 0,l,...,oo, then we may write w(x) as 



w(x) 



— J2 {l(^2|m 2 )| 2 (x v x b \e- s ^ + rt +2m+ ^ + ^\x v x b ) 

S m=0 



( e ^2 _ e -^ 2) _ \( x ^\ mv )\^ XtXb \ e -°iP 2 +pt+2 m +h\\]\ XtXb )( e 4v\\ _ e -«NI) 

(55) 
It should be obvious then that each term in w(x) in in fact a function of only 
one coordinate, x 2 . Indeed, one may even check that each term in the series 
above may be regarded as a 1 + 1 dimensional determinant times the square 
of the amplitude for the corresponding harmonic oscillator mode. 

The vacuum energy density in terms of the eigenstates of harmonic modes 
b\b 2 and c^c. To obtain the total vacuum energy, W, we have to integrate 

12 



over the phase space. After this step the sum over m can be performed and 
we obtain, 



L 2 L 2 

W = — limu> £ = — lim 

2ir e^o 2rr e^o 



J_^s(^+u b -E 2 ) coth ( SjB2 ) _ C oth(s||77| 



(56) 

The L 2 factor appears because we have considered the space volume to 
be a square box. This equation expresses the vacuum energy of a fermionic 
system with a domain wall that can be approximated linearly near the defect, 
and in the presence of a constant electromagnetic field, for a particular con- 
figuration (F ■ 77 = 0). We can see that in the case where the electromagnetic 
field is absent the vacuum energy vanishes, which is consistent with our non 
interacting fermions normalization condition. 

The regularized expression for W diverges as e goes to zero. The di- 
vergences in Eq.(|56|) should be isolated in some terms as a first step and 
afterwards suppressed by some renormalization procedure. To determine the 
divergent parts in e, we use the expansion, 






(57) 



valid for u 2 < tt 2 , B 2 k, being the Bernoulli numbers. We obtain, 



7-2 42 



r(-2,eA)-r(-2, 



-kA s 



+ B 2 E 2 



TT A 

r(o,eA)-r(o,— ) 

t, 2 



+ 



+ 



^ 2 2k B 2k E 2 2 k - 1 A 2 - 2k 



k=2 



(2k)\ 



TT A 

T{2k-2 1 eA)-T{2k-2 1 — ) 

E2 



\m 

2tt 2 



2|r/|e 2 



log 



71 



l\V\ 



-B 2 \V\ 



+ B 2 \r]\ log e — 



£ 

fc=2 



2 2k B 2k rr 2 



(2k)\ 



where 
F 1 (w v ,w b ,E 2 ) = 



(58) 



dss e 



-2-sA 



coth(s£ 2 ), F 2 (\t]\) = / dss' 1 coth(s|r/|) , 

(59) 
A = w v + Wb — E 2 and T(a,x) is the incomplete Gamma-function. In the 
expression (p8|) singularities are still present in the terms T(—2,eA) and 



13 



r(0, eA). To isolate them we use the formula, 

oo 

E^7> ( 60 ) 



oo j c 

T{a,x)=e- X x a * 



n=0 



where U^ are the Laguerre polynomials. We get at the limit e — > 0, T(— 2, eA) ■ 
2A 2 £ 2 and r(0, eA) — » J^ ^j-j-. This last series diverges logarithmically. Af- 
ter some rearrangements, we obtain, 



W e 



r2 1 11 oo 

z ' " 71=0 

^r(-2,^)-B 2 £; 2 r(o,^)-iog(^)-*N 

-C/2 &2 &2 \V\ 

.V A 5 2fe ,2i^ ^^ _ ^ eA) _ r(2fc _ 2) vr^ 



£7- 



^ y^ -P2fc /ffg\2fc 

E 2 ^(2fc)! 1 A J 



fc=2 



(61) 



The two first terms in flB"T|), contain quadratic and logarithmic divergen- 
cies which may be suppressed by a renormalization procedure. It may be, at 
first sight, shocking to see that the divergent part of W e is in fact not a finite 
degree polynomial in the external field F^ and its derivatives, as the usual 
divergence theorems for a standard Quantum Field Theory imply ||. The 
reason for this seemingly contradictory result is that the hypothesis for those 
theorems do not hold in the present case. Firstly, the mass term of the 'free' 
fermionic field is not a constant, thus the inhomogeneity of space already 
changes the main assumptions, like power counting behaviour. Note that in 
our case there is no unique large momentum behaviour for the propagator. 
And secondly, the would be external field F appears (after diagonalization) 
in such a way that it plays a similar role to the inhomogeneous mass. Then 
the external field is more likely to appear in a similar way to a mass in a 
standard divergent expression, which is indeed a non-polynomial dependence. 
In spite of this, a renormalization prescription may indeed be used for this 
quantity. For example, we may realize that W e , if regarded as a function of 
E2, has divergences where Ei appear as a pole, a constant or a linear term 

14 



in E 2 . Thus we may cancel all the divergences of this system by including 
three counterterms, namely, by adding to W e the 1-loop counterterm action 
SW e defined by: 

8W e = a" 1 {EiY 1 + a + a x E 2 (62) 

where the a/s are divergent constants. These three constants require the 
use of some renormalization conditions to fix them; in our case one could of 
course use the Laurent expansion of the actually measured W in order to 
fix those coefficients. It is amusing to note that this procedure requires the 
knowledge of the full dependence of W e on F, since the divergent part is not 
just a single polynomial in F. 

4 The general case 

Let us call Fn and F± respectively the projections of F onto the direction of 
7] and onto a direction orthogonal to rj (with respect to the scalar product 
defined in section 0) Then the square modulus of the Dirac operator for 
general configuration of the external vectors, can be written in the form, 



n = -d 2 +j{x 2 (Fl+F^)-[x.(F ± +F ll )} 2 }+(x.r]) 2 -ixxd.(F L +F ll )+-f.(F-r ] ). 

In the same way as before we define versors fj, t and b, where now £ M is a 
unit vector in the direction of F±, i. e. in the direction 

F M - F ■ f)rj»- (64) 

We can now define creation and destruction operators a a and a^, anal- 
ogously as it has been done in Eq. fli"5D , with a = r],t,b, and u> a defined 
by 






2 ± -L , 2 2 2 -L / C r\ 



In terms of these, H. can be written as, 



oo (uj 2 — UJ 2 ) 
H = ^uo a (2a) a a a + 1) - 4(^^ ^) 1/2 (a t a^ + a t a\ + a\a v + a\al) 



uJ„ 
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i(uj 2 b ■ 


-^ 2 ) 1/2 ( 
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-a t aj). 
(66) 



We see that 7i is again a bilinear form in the creation and destruction 
operators a a and a' a : 



H = a-Ajiai cijBjiai H — a'jB^aJ + constant 



2 -v-j.--. ■ 2 ~J i' 



(67) 



where Aji and .B^ are components of 3 x 3 matrices, A Hermitian and B 
symmetric: 




A 



with the constants defined by 



1Mb . ^t, 
h J— )u t s 

u t V Ub 
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(68) 
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UJt 



LOT 
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^6 - ^t 



(69) 



2? 

We also see that in this case it is possible to use a transformation, defined in 
terms of 3 x 3 matrices U and V, such that the Hamiltonian is diagonalized. 
The properties of this transformation can be summarized by the equations: 



U ] U - V*V = I U l V - V l U 



and 



AU + %B*V = UE 



iBU + A*V 







-VE . 



(70) 



(71) 



As in the simpler 2x2 case, one may obtain linear equations that deter- 
mine the transformation matrices U and V. In order for those equations to 
have a non-trivial solution, it is necessary to demand the condition: 

det [A*(B*)~ 1 A -B + {(B*)~ 1 A - A*{B*)~ 1 )E - {B*)~ 1 £ 2 ] = (72) 
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where E are the 'energies', i.e., the numbers appearing in the diagonal of the 
matrix E. Equation (|72| ) may be explicitly solved, what yields 6 solutions. 
Of these 6 solutions, we should eliminate 3 spurious ones, since they are not 
consistent with the algebraic relations that are verified by U and V. This 
leads to the following three solutions. One of them vanishes: E\ — 0, and 
the other two have a rather cumbersome expression, which may be written 
in terms of the previously defined u parameters as follows: 

S 2 = V2 {u 2 + u 2 + 3u 2 b - {u-^iul - ^ 2 )(64^ 3 + 96^^^)] 

W v - Uuy t + 13^ 4 + 6(c 2 - ^l + 9w£)])*}* (73) 

£ 3 = V2 {ul + u 2 t + 3u 2 b + (lu- 1 [{u 2 h - u 2 ) {U^ t + 96a; t 2 ^V^)] 

+^ - lAuy t + 13^ t 4 + 6(^ - co 2 U + 9^)])^ • (74) 

^From these solutions the Hamiltonian operator can be written in diag- 
onal form and the determinant of the Dirac operator can be obtained for 
general F and 77. Calculations in this case are more involved, but following 
along the same lines as we have done in section || it should be possible to 
obtain the vacuum energy in the general case. Progress on this subject will 
be reported elsewhere. 
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